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Using non-equilibrium Brownian dynamics computer simulations, we have investigated the steady 
state statistics of a polymer chain under three different shear environments: i) linear shear flow in 
the bulk (no walls) , ii) shear vorticity normal to the adsorbing wall, iii) shear gradient normal to the 
adsorbing wall. The statistical distribution of the chain end-to-end distance and its orientational 
angles are calculated within our monomer- resolved computer simulations. Over a wide range of 
shear rates, this distribution can be mapped onto a simple theoretical finite-extensible-nonlinear- 
elastic dumbbell model with fitted anisotropic effective spring constants. The tails of the angular 
distribution functions are consistent with scaling predictions borrowed from the bulk dumbbell 
model. Finally, the frequency of the characteristic periodic tumbling motion has been investigated 
by simulation as well and was found to be sublinear with the shear rate for the three set-ups, which 
extends earlier results done in experiments and simulations for free and tethered polymer molecules 
without adsorption. 

PACS numbers: 



1. INTRODUCTION 

Adsorption of macromolcculcs in flowing fluids arises 
in countless practical or potential applications. For in- 
stance the process of implanting bio-tissues happens in 
the shear flow field of the blood stream and represents a 
key challenge f]\ in the cure of disabilities due to malfunc- 
tioning organs. In gene therapy, the transport of DNA to 
a desired target site is an another big challenge. [1] The 
control of proteins changing their biofunctionality during 
adsorption to surfaces 3] is also a quite interesting topic 
for biotechnical application. Furthermore, the influence 
of adsorption under flow is important in microfluidic de- 
vices. [Jl For example, shear flow generates stretched 
polymer chains which can be used for microcircuit con- 
struction from single molecules, d, 0] 

Experimentally, polymers under shear flow have been 
investigated for many decades. 0, IE i, [M [H [H El, 
[3, [H, m, [13, Ell Due to both an elongational and a 
rotational components stemming from the flow inhomo- 
geneity, the dynamics or statistics of sheared molecules 
is complex. Tumbling motion was recently addressed ex- 
perimentally in the bulk, [H,!!!,!!!! but not for adsorbed 
polymers. 

On the theoretical side, Hinch has developed a the- 
ory for a continuous flexible string under shear flow, (2^ 
and also examined the effect of weak Brownian motion 
on the polymer. Recently, analysing the statistics of 
the chain end-to-end extension and orientation, Chertkov 
and co-workers have addressed the polymer tumbling mo- 
tion under shear flow. [24] Furthermore, Winkler has an- 
alytically studied the dynamics of semiflexible polymers 
under the influence of shear flow. 25] Numerical studies 
based on a simple polymer dumbbell in linear shear flows 
have been carried out to analyze the cyclic motion by 
Celani et al. [11] and Puliafito et al. 27]. As far as com- 



puter simulations on the tumbling motion are concerned, 
free chains [19, 28, 29] and end-tethered polymers (30l.[3l| 
and confined polymers in microchannels [s^ ] were consid- 
ered. 

Although the statistics of the polymer end-to-end ex- 
tension and angular orientation have been investigated 
experimentally and theoretically in the past, most of the 
relative work is essentially focused on the bulk case. In 
this paper, our main goal is to address these properties 
for an adsorbed polymer chain under linear shear flow. 
For a pulled adsorbed chain, a simulation study was re- 
cently performed in ref [33| . The problem of the ad- 
sorbed case is more complex than in the bulk since ad- 
sorption will interfere and compete with tumbling motion 
in the shear flow. Moreover, the chain conformational 
behaviour will depend on the direction of the shear flow 
relative to the adsorbing substrate. For a planar adsorb- 
ing wall, there are two basic possibilities for the shear 
flow direction: the shear gradient can be directed either 
perpendicular or parallel to the adsorbing wall. 

In this work, we consider these two cases in conjunc- 
tion with the bulk system as a reference and calculate the 
probability distribution function (PDF) of the chain ex- 
tension by using monomer-resolved Brownian dynamics 
computer simulations. We flnd that the PDF of the poly- 
mer chain extension can be mapped by a simple theoret- 
ical flnite-extensible-nonlinear-elastic (FENE) dumbbell 
model via two fit parameters, namely an effective shear 
rate and an effective spring constant. The tails of the 
angular PDFs obtained by the simulations are in good 
agreement with theoretical predictions stemming from 
the bulk. Surprisingly, they fit even better than in the 
bulk for a wide shear rate regime when the shear gradient 
is directed parallel to the adsorbing wall. For the situ- 
ation where the shear gradient is perpendicular to the 
wall, disturbances of the bulk fluctuations due to wall 
adsorption lead to significant deviations from a simple 
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dumbbell model. We also explore the inner bond-bond 
angles of the chain and the frequency of the character- 
istic periodic tumbling motion which was found to be 
sublinear with the shear rate. This extends earlier ex- 
perimental and simulation studies for free and tethered 
polymer molecules without adsorption. 

The outline of this article is as follows: Section [2] de- 
scribes the details of the simulation model. Section [3] 
contains the results, in which section [3T] deals with the 
statistics of the chain extension; 13.21 shows the angu- 
lar degrees of freedom (for in-adsorption-plane and off- 
adsorption-plane orientations, respectively); 13.31 is dis- 
cussing the shear rate dependency of the bond-bond an- 
gle; and the last part 13.41 is devoted to analysis of the 
tumbling motion frequency. Finally, a summary of our 
work is presented and discussed in section U] 



THE MODEL 
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FIG. 1: Sketches of simulation model: (a) for the shear gra- 
dient vorticity normal to the wall (NV); The wall normal is 
along the direction, (b) for the shear gradient normal to 
the wall (NG); (c) for the bulk (linear shear flow without 
walls). Also shown in (c) and (d) are the angles 9 {6') and 
{(j)') that characterize the end to end vector. 

This section provides details about our polymer chain 
model employed in non-equilibrium Brownian dynamics 
(BD) computer simulations. [H, [IHl The macromolec- 
ular chain is made up of a linear sequence of A'^ = 48 
coarse-grained bead-spring monomers, which are labeled 
i — 0, 1, 2, • • • , iV — 1. Thereby, the equation of mo- 
tion of the i-th monomer located at a position ri{t) = 
{xi{t),yi(t), Zi{t)) at time t is given by the finite integra- 
tion scheme 

r,{t + 5t) = Y,{t) + -^F.St + 6G^ + jq^it)Ste,, (1) 
kB-l 



where {t + 6t) is the updated bead position at a later 
time t + St. Hydrodynamic interactions are neglected. 
[3^ We now explain in details all the constitutive terms 
of Eq. ID). Dq, ks and T stand for the free diffusion 
constant, Boltzmann's constant and the absolute tem- 
perature, respectively. The vector is the total con- 
servative force that will be fully described later. 5Gi is 
the Gaussian stochastic displacement with zero mean and 
variance 2Do6t for each Cartesian component. The last 
term in Eq. [T] stems from the contribution of the shear 
flow being in the a;-direction, with e^, denoting the corre- 
sponding unit vector, and 7 the shear rate. qi(t) in Eq. [1] 
can be either yi{t) or Zi{t) depending on the chosen shear 
gradient direction. More precisely, (i) qi{t) = yi{t) refers 
to a normal vorticity (NV) direction perpendicular to the 
adsorbing wall, as sketched in Fig. [TJa); (ii) qi{t) = Zi{t) 
corresponds to a normal gradient (NG), see Fig. [IKb). 
In other words, the shear gradient is in-plane with the 
adsorbing wall for the NV case whereas it becomes off- 
plane for the NG case. The aforementioned conservative 
force Yi has three contributions: 

1. Steric effects (monomer- monomer and wall- 
monomer repulsive interactions) are taken into 
account via a truncated and shifted purely 
repulsive Lennard- Jones potential of the form 



ULj{.r) = 4e 



r J \ r / \rc J \rc 

(2) 

with a cutoff at its minimum Tc — Here 
the bead diameter ah and the strength of the in- 
teraction e = ksT represent the energy and length 
units, respectively. 

2. The "spring" is modeled via a finite extensible 
nonlinear elastic (FENE) [l^l potential to ensure 
the connectivity between adjacent beads along the 
backbone, which is given by 



Up EN sir) 



Rq 



(3) 



where the FENE cut-off (i.e. the allowed maximum 
bond length) is set to Rq — l.Scrf, and the spring 
constant is if = 27e/ab^ . 

3. To mimic the adsorption process, a strongly attrac- 
tive wall defining the (x, y)-plane at z = is consid- 
ered. Thereby, we employed a Van der Waals-like 
attractive potential 



Uads{z) 



-A 



(4) 



with ^0 = 5- A more detailed modeling for a real- 
istic surface was used e.g. in ref [38| . 

The BD time step is set to St = 2 x 10~^r, where 
T = (jI/Dq sets the time unit. Typically 5 x 10^ BD 
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time steps are used for data production, during which the 
first 5 X f 0^ steps are for the chain relaxation. 1^°^^^ i) 
Simulation snapshots are viewed in Fig. [2] 



3. RESULTS 
3.1. Statistics of the polymer extension 

We consider a single polymer chain which is advected 
by a shear flow and stretched by velocity inhomogeneity. 
The degree of the polymer stretching is characterized by 
the chain's end-to-end distance Re, for which we have 
calculated the probability distribution function (PDF) 
under different shear flow geometries. 

The normalized PDF of the rescaled size extension 
Re/Rm, where Rm = {N — l)Ro is the fully stretched 
chain length, is plotted in Fig. [3] for three typical se- 
lected shear rates (-yr = l(n), 10(o), 50(A)) for three ge- 
ometries: (a) bulk, (b) NV, in-plane shear gradient and 
(c) NG, vertical shear gradient. 1'^°*'' 2] 

A general feature is that upon increasing the shear 
rate, chain stretching is always favored. In terms of size 
distribution, this shear-rate induced stretching manifests 
itself either as a shift of the peak to larger sizes (see Fig. 
1311a, b)), and/or as a thickening of the tail (see Figl^Kc)). 
The strongest stretching, at prescribed reduced shear 
rate 7r, is obtained for the NV case. This feature can 
be simply explained as follows. The NV case, due to 
the rather strong adsorption, leads to a chain swelling 
(in the lateral dimensions) compared to the bulk confor- 
mations already at equilibrium (i.e., without shear flow). 
To elucidate this idea, we have also systematically shown 
the purely two-dimensional case (filled symbols) in Fig. 
|3l^a,b) as a reference. Moreover, for the NV case, the 
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FIG. 2: Typical snapshots of polymer chains under linear 
shear flow in a;— direction. The shear gradient direction is nor- 
mal to the wall-plane (NG). Here the three columns (a)(b)(c) 
show configurations for three shear rates 7t = 1, 10, 50. 



[1] Given the strong adsorption, the NG case (normal shear gra- 
dient) requires more computational time for good statistics, so 
that even 10® BD steps were used. 

[2] In the figures of the results section, all the error bars were cal- 
culated using five (for the bulk case) and ten (for NV and NG) 
independent simulation runs. 






□ 1 

Yi= 10 
A -p: = 50 

— K_= 1.5K= 15E/a^",yi = 7 

— K_= 1.5K= 15E/a^",-iT = ™ 

— K^= 1.5K= 15e/o^", 17 = 350 




■ 

NV 
















0.2 


0.4 0.6 %7 





R/R 



R/R 





□ YI=1 

■ Yi=iO 


NG 




A yr = 50 -, 
■ K^ = l.SK=15E/a^".yt = 7 
K^= 1.5K= lSE/a^-,Yi = 70 
K^ = l.SK=15E/a'.-jT = 350 

— K = 5K = IQOda^-, yr = 7 

— K =5K = 200E/a^l-iT = 70 

— K =5K = 2Q0E/a-.Yr=350 


if 

0.2 





R/R 



FIG. 3: The probability distribution functions P{Re/ Rm) of 
the reduced end-to-end distance Re/Rm for (a) the bulk, (b) 
the normal vorticity (NV) case and (c) the normal gradient 
(NG). The symbols □, o, A stand for the shear rate 7r — 
1,10,50, respectively, where in (a) and (b), the pure two- 
dimensional (2D) case is shown as a reference (filled symbols). 
For the FENE-dumbbeU model (curves), the chosen fitting 
spring constant(s) parameter(s) is(are) indicated. 



shear flow deforms the polymer chain much more effi- 
ciently than in the NG case, where the shear rate gradi- 
ent direction is off plane. Clearly, deformations occur in 
the NG case because of temperature-induced fluctuations 
in the heights of the monomers. 

To broaden our understanding of the statistics of 
the polymer chain under shear flow, we have addition- 
ally considered the well known so-called FENE-dumbbeU 
model. The equation, describing the evolution of the 
end-to-end vector Re of a FENE-dumbbeU (i.e., a dimer), 
is given by 



Re(t + (5t) = Re(t) 



Do 



FpENESt + 6G + jQ{t)Ste^ 



(5) 

where Fpene — -^UpENEiRe) is the FENE force; 5G 
is the Gaussian term as in eq. [l] Q is the y— (or z—) 
component Rey (or Rez) of Rg assuming the shear gra- 
dient applied in the y— (or z~) direction. To take ad- 
sorption into account in the most simple and physically 
sound manner within the dumbbell model, we have in- 
troduced anisotropic FENE-spring constants such that 
= Ky ^ Kz ^ K in the bulk and [K^ > {K^ = Ky = 
K)] for NV and NG cases. The stochastic equation ([5]) 
was solved by using a straightforward one-particle BD 
scheme integration. One has to bear in mind that this 
dumbbell approach is certainly relevant in the limit of 
strong stretching, where the degrees of freedom of "in- 
ner" monomers of a real polymer chain would be indeed 
heavily reduced, but can not be suitable at low shear 
rates where the real chain size is dictated by the compli- 
cated distribution of the monomers. 
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FIG. 4: The first three moments of the probability distribu- 
tion function of end-to-end distance (see Fig. |3} for the three 
cases (NG, NV and bulk) plotted against the shear rate: (a) 
the mean value of the end-to-end distance (|-Re|); (b) the vari- 
ance (cr^); (c) the third standardized moment (S), so-called 
skewness. Here the solid symbols stand for the purely 2D 
results. 



By choosing a suitable choice of parameters (effective 
spring constant (s), effective shear rate) for the dumb- 
bell model, one could always perfectly match any PDF 
of the fully monomer resolved chain (results not shown) . 
However, in order to make the comparison more physi- 
cal between the dumbbell model and the fully monomer 
resolved chain, we kept the spring dumbbell constant 
fixed as well as the ratio between the shear rates, see 
Fig. El (7 : 70 : 350 = 1 : 10 : 50). Given the very 
crude representation of the chain via a dumbbell, Fig. 
[3Ua,b) shows a remarkable good agreement between BD 
data and the dumbbell model, especially at high shear 
rates and/or large chain size. The good agreement be- 
tween these two approaches in Fig. ^h) for the NV 
case also demonstrates that the idea of mimicking the 
adsorption with a larger dumbbell-spring constant in the 
z-direction is fruitful. Nonetheless, the NG case is too 
complex to be captured by our simple dumbbell model. 
Indeed, physically, the deformation in this situation is 
due to the z-fluctuations of the chain-monomers which 
can not be suitably approximated by a dumbbell. How- 
ever, if we consider a stiffer (effective) dumbbell with 
Kz = 5K — 200 for the NG case, we can recover again a 
good agreement with the BD data, see Fig. ^c). 

To further characterize the statistics of the adsorbed 
chain under shear flow, we have calculated the flrst three 
moments of the end-to-end size distributions (the PDF in 
Figure[3|), which are (i) the mean value (|Re|) = (-Re), (ii) 
the variance ct^, and (iii) the normalized third moment 
(so-called skewness, S). 



The mean value of the end-to-end distance is obtained 
as follows 



(|Re|> - {Re) 




E 



M 



^1 Ir^v-i - rol 



M 



(6) 

where Rem is the absolute value of the end-to-end vector 
in the m— th configuration and the sum runs over the to- 
tal number of the configurations M. Figure Ufa) shows 
that (Re) grows with increasing (7T), as expected. Nev- 
ertheless, in the NG case, the shear rate dependency of 
the end to end distance is relatively weak compared to 
the other flow geometries. 

The second moment, which describes the degree of 
broadening about the mean value, is defined as follows 



Em=l i^e 



{Re)r 



M 



(7) 



Figure IHb) clearly shows that the variance grows 
with increasing shear rate. This effect is actually due to 
the tumbling (more on that later) allowing the chain to 
cyclically stretch and recover its coil-like conformation. 
Thereby, the higher the shear rate the more stretched the 
chain can be, and consequently the larger the variance 
becomes. In the NG case, where stretching at prescribed 
shear rate is much weaker, the broadening of the size 
distribution is reduced accordingly. 

The skewness S, a central quantity in probability the- 
ory and statistics, [s^ is defined as S = fis/cr^, where 
fi is the third moment with respect to the mean, and a 
is the standard deviation. Explicitly, the skewness was 
computed as follows 



rr3 



M Em=l (-^e 



{Re)f 



M Em=l {Reva {Re)) 



3/2 



(8) 



It provides a quantitative and simple measure of the 
asymmetry of the probability distribution function. 
Hence, the qualitative difference between the PDF in the 
NG model (see Fig. EJc)) and those found in the other 
cases (see Fig. [3l^a)(b)) can be rationalized in terms of the 
value as well as the sign of S, see Fig. Hfc). For the shear 
gradient normal to the wall (NG), Fig. [DJc) shows that S 
is positive and increases with shear rate. The positivity 
of {S) is a signature of the statistical preponderance of 
the coil-sized chain configurations in the NG model. On 
the other hand, the skewness becomes negative and de- 
creases (i.e., its magnitude increases) with growing shear 
rate in the NV model. This feature clearly indicates the 
preponderance of stretched chain conformations. 



3.2. Statistics of the polymer end-to-end vector 
orientations 

We now address the angular degrees of freedom of the 
end-to-end polymer chain vector, namely (j) and 6 (de- 
picted in Fig. [He)), or (/)' and 0' (depicted in Fig. md)). 
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FIG. 5: The PDF of the end-to-end vector angles (see Fig. 
[ijc)) for the three cases (bulk, NV and NG). The symbols 
and the curves stand for the same parameters as those usee 
in Fig. [31 



with respect to the prescribed shear flow x-direction a: 
well as the shear gradient direction. Whereas 4)/(j)' corre 
sponds to the usual azimuthal angle, 9/9' is not the pola 
angle but rather the angle formed between i?e and th( 
plane containing the shear- and the gradient-direction 
see Fig. [llc,d). Thereby we have kept the same (well es- 
tablished) definitions used in earlier publications. [23.l40| 
The results for the angular distribution functions are 
shown in Fig. [5l As far as the azimuthal angle (f) dis- 
tribution is concerned, BD simulation data as well as 
the dumbbell model show a strong probability peaks in 
the vicinity oi ((> — Q, tt, 27r in the bulk (Fig. [S^a)) and 
the NV case (Fig. Efc)). Interestingly the corresponding 
peak heights are not identical and, in fact, decrease with 
growing angles. The simple underlying physical mecha- 
nism accounting for this effect is that the torque, (result- 
ing from the shear gradient [T9|). exerted on the chain 
favors more the conformations with cj) Q than those 
with (j) — > 27r. This finding of the asymmetry in Pi^cj)), al- 
though unexplained, was also reported in the earlier work 
of Celani et al. [2^. Moreover, the larger the shear rate 
7T, the narrower the density probability becomes around 
the flow (x-axis) direction, see Fig. [3a,c). However, the 
PDFs are somewhat different for the NG case, see Fig. 
[5]Je), due to the shear gradient direction being off-plane 
hereby. 

As far as the 0— PDFs are concerned, the shear flow 
induces a concentrated distribution around 9 = for all 



FIG. 6: The PDF of the in-shear-gradient-plane angle (p {(j)') 
(sketched in Fig. [TJc) and in Fig. [Hd)). The solid curves 
are the theoretical scaliiig predictions sin^^0 or sin^^<jf>') 
of the dumbbeU model. [M [H 
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FIG. 7: The PDF of the off-shear-gradient-plane angle {8') 
(sketched in Fig. (TJc) and in Fig. H^d)). The solid curves 
are the theoretical scaling predictions or d'~^) of the 

dumbbeU model. [H,!!! 



shear flow geometries, see Fig. [S][t>,d,f). 1'^°*'' 3] -pj^g 
narrowness of the peak increases again with the shear 



[3] Note that the symmetry of the shear flow imposes a symmetry 
of P{e) about 9 = 0. 
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rate, as expected. 

In a general manner, the agreement between our simu- 
lation data and the dumbbell model (see Fig. [5]) is not as 
satisfactory as that found for the size distribution (com- 
pare with Fig. [3|). This is quite acceptable, given the 
crudeness of the dumbbell model. Nonetheless, the es- 
sential feature of strong depletion zones in the angular 
PDFs, are qualitatively well captured by this dumbbell 
model. 

The angular distribution functions were also recently 
studied analytically by Chertkov et al. [24[ for a linear 
dumbbell in the bulk. It was thereby found to obey the 
following scaling law [l^l P{<P) ^ sin~^ 0, which was ex- 
perimentally confirmed by Gerashchenko and Steinberg. 
[2l| To establish a comparison between our simulation 
data and this theoretical prediction, the PDFs are plotted 
for a wider range of the shear rate on a linear-logarithmic 
plane in Fig. [6] A good agreement is found between 
the simulation and the scaling prediction. This is a re- 
markable result, since a simple harmonic dumbbell model 
can capture the correct scaling behavior of P{4>) of a fuF 
chain. 

Concerning the 9 distribution function, Chertkov e 
al. [l^l have also shown a scaling behavior for the tai 3,, 
of the form P{0) ^ within the framework of a lin 
ear dumbbell in the bulk. Comparison in the past witl 
experimental results [2lj demonstrated already an excel 
lent agreement. Our results for the 6'-PDF can be founc 
in Fig. [71 The scaling prediction is again in remarkabL 
good agreement with our BD data. Interestingly, the 
agreement seems to be even better for the NV case, where 
not only the tail is dictated by but the full (?-range. 
This suggests that the algebraic decay P{0) ^ 9^^ is even 
more relevant for (nearly) two-dimensional systems. 



3.3. Bond-bond angle statistics 




FIG. 8: The average bond-bond angle (a) plotted against the 
reduced shear rate 7T. Inset: The order parameter P2(cosa). 
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To characterize the intra chain property under shear 
flow, we have considered the bond-bond angles (BBAs). 
This is a complementary approach to investigate the 
chain stretching and coiling features. The results are 
depicted in Fig. [51 In a first low shear rate window 
('7''' ^ 20), there is a rather strong shear rate depen- 
dency for the NV, bulk, and 2D cases. Typically, this de- 
scribes the relevance of coiled conformations. At higher 
shear rates, this dependency becomes much weaker and a 
plateau-like regime (about tt) is observed. This feature, is 
the signature of strongly stretched polymer chains. As far 
as the NG case is concerned, the shear rate dependency 
is much weaker as expected. This result is fully consis- 
tent with the positive skewness reported in Fig. ^c) for 
a wide range of the shear rate. A complementar y ob serv- 
able provided by the following order parameter [4l| P2 is 
also shown as an inset in Fig. [51 



FIG. 9: The mean local bond-bond angle (ai) of monomer 
(j). The solid symbols in (a) and (b) represent the purely 2D 
case. 



P2 (cos a) ~ 




We now would like to address the local bond-bond an- 
gles, which physically characterizes the degree of local 
chain stretching. The mean BBA (ai) is plotted vs. the 
monomer position index (i) in Fig. [51 As expected, a 
symmetry around the chain center is revealed, see Fig. 
[HI The two chain ends are always looser and the local 
stretching, is smoothly increasing upon approaching the 
center of the chain. Besides, the local chain stretching 
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FIG. 10: Reduced tumbling frequency vt as a function of the 
reduced shear rate 7T. The lines correspond to the best fit. 

increases with growing shear rate, see Fig. |21 



3.4. Tumbling 

The tumbling process of the chain under shear flow 
is the mechanism by which energy is transferring, or in 
other words, the polymer cyclically absorbs (stretching) 
and releases (collapsing) energy from and to the sur- 
rounding fluid. This particular polymer motion in shear 
flow has been thoroughly studied in the past. p^.[20l.[30j 
All these studies address the bulk case or grafted 

polymer chains, [lOl but without considering the problem 
of adsorption. In this work, the frequency of this cyclic 
motion is calculated for the three geometries (bulk, NV 
and NG) via the Fourier transform of the time-dependent 
radius of gyration, which is 

Rg{t) = ^ ^^0'(l-»(^)^-I-C.„..(t))2 ^ ^^^^ 

where rc.m.(^) = T^f^o^ 1-^ is the chain center of 
mass at time t. The first main peak position of the re- 
sulting Fourier transform of Rg{t) (not shown here) is 
chosen as the tumbling frequency. 

Figure [TO] displays the reduced tumbling frequency 
{vt) as a function of the reduced shear rate 7r. The 
purely 2D case is also shown as a reference. Thereby, 
scaling behaviors of the form 

are observed , with A = 0.90,0.71,0.72 for NG, NV, and 
bulk, respectively. Hence, the bulk and NV cases, exhibit 
(within statistical uncertainties) an identical scaling be- 
havior for the tumbling frequency. The value found for 
the exponent A « 0.7 is in very good agreement with the 



one reported by Schroeder et al. [T^ (A = 0.67), based on 
the angular power spectral frequency, for Brownian rods 
in the bulk. As expected the NG case exhibits a different 
scaling behavior (here A = 0.9) due to the chain adsorp- 
tion that weakens the effect of the velocity gradient in 
the normal direction. 

To understand the physical origin of the observed scal- 
ing behaviour of the tumbling frequency in Fig. 1101 we 
are goiiig to employ the arguments advocated by Teixeira 
et al. [20|. By decomposing the tumbling motion into 
typical time periods: (i) a stretching phase {tstretch), (h) 
an alignment along the shear direction {t align) i (hi) a flip 
motion {tfup), and (iv) a collapse event {tr.niianxe), we can 
write an expression for the tumbling frequency [l3 . [20| v 

T 

VT OC . (11) 

^stretch ^a/zgn "t" tflip ~\~ tcollapse 

The three processes of polymer stretching, aligning and 
collapsing being advection-driven, Teixeira et al. (20| 
have shown that tstretch, taUgn, and t collapse scale like 
(7T)(fe) ' ^^cre {X) is the average chain extension in the 
x-flow-direction and {62) is the chain extension in the 
shear gradient direction. On the other hand, the flipping 
process is a diffusive motion such that tfUp ~ ofs^) ^ 
(62)^^^, [lOl where D{S2) is the diffusivity in the shear 
gradient direction and found to be D{62) ~ S2 "^^^ from 
BD simulations by Doyle et al.. Following the same 
procedure with our simulation data, we were then able 
to extract the time duration scalings which are gathered 
in table m 



TABLE I: Characteristic times involved in the tumbling mo- 
tion discussed in the text. 





tstretch 


talign 


tcollapse 


t flip 


bulk 

NV 

NG 


~ (7r) 

- (jr)-'-'' 

- iir)-"-'' 


~ (7r) 

- (7r)-«-«* 

~ (7r)-«-«^ 


~ (7t) 

~ (7r)-»-«^ 

- iir)-"''' 


- (7r)-»'" 
- (7r)«'« 



These results in Table U are in good agreement, for 
the NV case and the bulk, with the observed exponent 
(A = 0.7, see Fig. [TOl). Even for the NG case, the scal- 
ing predictions are validated for the convective processes 
(stretching, alignment, collapsing). Nonetheless, it fails 
for tfiip in the NG case because the strong adsorption 
weakens the diffusivity in the shear gradient {z—) direc- 
tion D{62). 

4. CONCLUSION 

Using Brownian dynamics simulations, the statistical 
properties of adsorbed polymers under a linear shear flow 
have been studied for the two cases: in- adsorption-plane 
shear gradient (or equivalently the shear vorticity nor- 
mal to the adsorbing wall), and the shear gradient per- 
pendicular to the wall. The three-dimensional bulk and 



8 



the purely two-dimensional cases under linear shear flow 
have been considered as well as reference cases. 

We have compared the behaviour of the full polymer 
chain to that of a simple dumbbell model. A good agree- 
ment is reached for the PDF of the chain end-to-end dis- 
tance. Besides, the theoretically predicted scaling be- 
haviour [U, [23 of the angular PDF tails is confirmed 
also in the case where the shear gradient is parallel to 
the walls (NV case) . In the opposite case of perpendicu- 
lar shear gradient (NG case), agreement is only found at 
high shear rates. It is remarkable that a simple dumbbell 
model can capture the correct scaling behaviour of the 
angular PDFs' tails of a full monomer-resolved chain. 

The chain coiling and stretching degree was character- 
ized by the inner bond-bond angles. For the NV, bulk 
and 2D cases, the shear rate dependency is found to be 
rather strong in the first low shear rate regime, but is 
much weaker at higher shear rates, and a plateau-like be- 
havior is approached once the shear rate is high enough. 
For the NG case, this dependency is much weaker for a 
wide range of shear rates. 

The tumbling frequency in the bulk scales sublinearly 
with shear flow rate, which is in agreement with the early 
studies. [H, [131 We found this sublinear scaling also in 
the adsorbed cases (NV, NG). The theoretical scaling 
analysis based on the dumbbell model, was confirmed 



for the NV case but not for the NG case. 

As an outlook, understanding the stretching and coil- 
ing mechanisms of polymer in flow leads to an im- 
proved control over chain conflgurations in various cir- 
cumstances of confinement, e.g. in microfiuidic devices. 
[3] Stretched polymer chains may provide valuable candi- 
dates for nanowires in microcircuits. Q Therefore a min- 
imalistic model, as the simple FENE dumbbell model, is 
useful in order to map the behaviour of a complex sys- 
tem and to predict the qualitative trends in various flow 
situations. 

For future studies, many polymer chains, i.e. a concen- 
trated polymer solution, may exhibit complex dynamical 
behaviour due to significant slowing down caused by mu- 
tual chain entanglements. |43|| It is further interesting to 
mix nanoparticles and polymers [4^ and expose them to 
different flow flelds and study their collective dynamical 
behaviour. [45] Finally, the influence of hydrodynamic 
interactions need more consideration, following the re- 
cent simulation method proposed in ref [36| . 
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